Astronomy & Astrophysics manuscript no. aa 


© ESO 2009 


November 25, 2009 





Determination of tlie Local Dark Matter Density in our Galaxy 



M. Weber' and W. de Boer' 



0^ 
O 
O 
(N 

> 
O 

in 

(N 

o 

u 

o 



Institut fiir Experimentelle Kernphysik, Karlsruher Insitut fiir Technologie (KIT), P.O. Box 6980, 76128 Karlsruhe, Germany 
Received 30 September 2009 / Accepted 9 October 2009 

ABSTRACT 

Context. The rotation curve, the total mass and the gravitational potential of the Galaxy are sensitive measurements of the dark matter 
halo profile. 

Aims. In this publication cuspy and cored DM halo profiles are analysed with respect to recent astronomical constraints in order to 
constrain the shape of the Galactic DM halo and the local DM density. 

Methods. All Galactic density components (luminous matter and DM) are parametrized. Then the total density distribution is con- 
strained by astronomical observations: 1) the total mass of the Galaxy, 2) the total matter density at the position of the Sun, 3) the 
surface density of the visible matter, 4) the surface density of the total matter in the vicinity of the Sun, 5) the rotation speed of the Sun 
and 6) the shape of the velocity distribution within and above the Galactic disc. The mass model of the Galaxy is mainly constrained 
by the local matter density (Oort limit), the rotation speed of the Sun and the total mass of the Galaxy from tracer stars in the halo. 
Results. It is shown from a statistical fit to all data that the local DM density is strongly positively (negatively) correlated with the 
scale length of the DM halo (baryonic disc). Since these scale lengths are poorly constrained the local DM density can vary from 0.2 
to 0.4 GeV cm"' (0.005 - 0.01 Mq pc"') for a spherical DM halo profile and allowing total Galaxy masses up to 2 • 10'^ Mq. For oblate 
DM halos and dark matter discs, as predicted in recent N-body simulations, the local DM density can be increased significantly. 

Key words. — Milky Way: halo, structure, dark matter, rotation curve, dwarf galaxies, dark matter mass, baryonic mass — 
Cosmology: dark matter, dwarf galaxies, structure formation 
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1. Introduction 

The best evidence for Dark Matter (DM) in galaxies is usually 
provided by rotation curves, which do not fall off fast enough 
at large distances from the centre. This can be understood, either 
by assuming Newton's law of gravitation is not valid at large dis- 
tances, the so-called MOND (Modified Newtonian Dynamics) 
theory (Bekenstein and Milgrom 1984; Bienavme et al. 2009), 
or the visible mass distribution is augmented by invisible mass, 
i.e. DM. For a flat rotation curve the DM density has to fall off 
like 1/r^ at large distances. 

Given the overwhelming evidence for DM on all scales from 
the flatness of the universe combined with gravitational lensing 
and structure formation we assume DM exists and try to con- 
strain the DM density from dynamical constraints, for which 
better data became available in recent years: 

1. the total mass of the Galaxy has to be about 10'^ solar masses 
dWilkinson and Evanslll999HBattagha et al.ll2005l: IXue et all 



120081) : 

the total mass inside the solar orbit is constrained by the 
well-known rotation speed of the solar system; 
the total matter density at the position of the Sun from 
the gravitational potential determined from the movements 
of local stars as measured with the Hipparcos satellite 
(Holmberg and Flynn 2004); 

the s urface density of th e visible matter at the position of the 
Sun dNaab and Ostriker|[2006. and references therein); 
the s urface density of the total matter at the position of the 
Sun (Kuiiken and Gilm orell99TllHolmberg and Flynnl2004l: 
[Bienayme et al. 2005); 



6. the shape of th e rotation curve within Galactic disc 
dSofue et"aLll2008l and references therein); 

7. the velocit y distr ibution above the Galactic disc (z > 4 kpc) 



cit y I 

(IXue et alJl2008h 

Towards the Galactic centre (GC) N-body simulations of struc- 
ture formation predict a DM density, which diverges at the cen- 
tre, thus developing a so-called cusp, with the slope varying 
more li ke 1/r than l/r^ dNavarro et alJll997t iMooreet al.|[l999t 
lRicottill2003h . On the contrary, rotation curves of nearby dwarf 
galaxies do not show such a cusp, but rather a constant density 
near the center, a so-called core (Oh et al J 120081; iGentile et al.l 
2007; Salucci et al. 2007). N-body simulations including the 
complicated baryonic disc formation indicated that the cusps 
may have been destroyed by the fluctuation in the gravitational 
potential of the baryons, which is muc h stronger than the poten - 
tial caused by the DM in the centre (Mas hchenko et alll2006l) . 
In addition, the DM may form a disc like structure by the in- 
fall of DM dwarf ga laxies, as shown b y recent high resolution 
N-body simulations dPurcell et alj|2009 ). Evidence for substruc- 
ture of D M inside the disc is supported by the structure in the gas 
flaring ( Kalbe rla et alj|2007l) . by the peculiar change of slope for 
the rotation curve inside the disc, which is not observed for stars 
outside the disc (Sect. 13.21) and by the structure in the diffuse 
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gamma radiation ("de B oer et alj|2005l) . 

Such substructure will not be investigated in this paper, but 
smooth DM halos with different (cored and cuspy) profiles will 
be compared with all available data. These will result in lower 
limits on the local DM density, since dark matter discs or other 
local substructure will only enhance the local density. 

A reliable determination of the local DM density is of great 
interest for direct DM search experiments, where elastic colli- 
sions between WIMPs and the target material of the detector are 
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searched for. This signal is proportional to the local d ensity. A 
review on direct searches can be found in the paper bv lSpoorien 
(I2007h . 

The structure of the paper is as follows: in section |2] the 
parametrization of the luminous matter and five different DM 
halo profiles are given. In Sect. [3] the experimental data used to 
determine the mass model are discussed. Then the numerical de- 
termination of the mass model parameters using a fit of all 
dynamical constraints is discussed in Sect. |4] At the end a sum- 
mary of the results is given. 

2. Parametrization of the Density Distributions 

In order to constrain the mass model of the Galaxy by data it is 
convenient to have a parametrization for both, the visible an dark 
matter density. These parametrizations are introduced here. 

2. 1. Parametrization of the Luminous Matter Density 

The density distribution of the luminous matter of a spiral galaxy 
is split into two parts, the Galactic disc and the Galactic bulge. 
The parametrization of the den sity distribution of the bu lge is 
adapted from the publication by lCardone and Serend (120051) 



Pb(r, z) = Pfc ■ — 



(1) 



For a good description of the RC near the GC the parameters 
of the bulge profile are found to be jb - /3h = 1-6, qb - 0.61, 
rt = 0.6 kpc, ro,b = 0.7 kpc and pb = 360 GeV cm~^, which 
corresponds to approximately 9.5 Mq pc"^. This high density 
might be influenced by the presence of black holes. At least one 
black hole w ith a mass of about 4.0 ■ 10^ Mq has been observed 
in the GC bv lReidI (|2009l) . 

The stellar contribution of the Galactic disc is split into two discs 
- a thin and a thick disc - which are usually parametrized by an 
exponentially decreasing density distribution. 

The para metrization of t he Galactic disc is taken from the 
publication bv lSparkd (l2007h 



Pd{r,z) = Pd ■ exp(-r/rd) ■ exp(-z/Zd)- 



(2) 



The parameter pd describes the density of the Galactic disc at 
the GC while id and Zd describe the scale parameter in radial and 
vertical direction. 

There is some freedom in the choice of the parameters for 
the Galactic disc. Its density in the GC is, as in case of the bulge, 
unknown, so it has to b e a free p arameter. For the scale radius 
we adopt the value from lHammer et al., (,2007i) 



rd = 2.3 + 0.6 kpc. 



(3) 



The scale height Zd varies for the different components: young 
stars are born near z = 0, so they have a much smaller scale 
height, while the older stars have been kicked around and reach 
scale heights up to 700 pc. The disc has two distinct populations 
with two different scale heights. The thin disc consists of young, 
bright and metal -rich stars which provide abo ut 98% of the to- 
tal stellar population jGilmore and Reidlll983l) . Its scale height 
was determined from star counts of stars with an absolute mag- 
nitude of My > 6 to be Zd ~ 270 pc (Kroupa et al. 1993). From 
measurements of the spatial distribution of K dwarf stars the ver- 
tical scale height turns out to be smaller than 250 pc (iRobin et~an 



Il996h . In the pubhcation bv lOiha etal] (Il996h the scale height 
of the thin disc was determined for stars with My > 3.5 to be 
260 ± 50 pc. Furthermore the interstellar gas and dust contri- 
bution is also located in the thin disc. The thick disc consists 
of old, metal-poor stars and could be the result of either an ear- 
lier thin disc or th e merging of a satellite galaxy with the MW 
( iRobin etalj|1996l) . Its scale height Ues between 700 and 1500 
pc and its local density is about 5% of the local stellar density of 
the thin disc. For that reason in this analysis the luminous matter 
is given only by the thin disc while small density contribution 
of the thick disc is neglected. For gas, which makes up 10% of 
the mass of the disc, the scale height varies with Galactocentric 
radius because of the decreasing gravitational potential at larger 
radii. Fortunately, the mass model of the Galaxy is not very sen- 
sitive to the exact value of the scale heigh t, so it is fixed at 32 
pc, which is close to the value adopted bv' Freudenreichl (Il998h . 
The mass of the bulge was taken into account, but the detailed 
structure of the mas s distribution (hole in the centre of the disc 
iFreudenreichI d 1 9981) . the bar structure and the black hole) are not 
taken into account in the parametrization, since the parameters 
of interest, i.e. the DM halo parameters, are insensitive to the 
detailed mass distribution in the centre of the Galaxy, since the 
gravitational potential in the centre is completely dominated by 
the baryonic matter. 

The parametrization of the visible mass discussed above 
leads to a mass of the Galactic bulge of about 1.1 ■ 10'" Mq. 
The mass of the Galactic disc varies for different fits because of 
the variation of the parameters pd and rd. It is in the range of 5 ■ 
10'° to 7 ■ 10'" solar masses. In addition to the luminous matter 
the density profile of the DM halo has to be parametrized. This 
is discussed in the following section. 

2.2. Parameterization of tine Dark Matter density 

The first a nalytical ana lysis of structure formation in the 
Universe by i GunnI (1 19771) predicted that DM in Galactic haloes 
are distributed according to a simple power law distribution 
p(r) oc r"''. However, later studies based on numerical N-body 
simulations (Navarro et al. 1997; Moore et al. 1999) found that 
the slope of the density distribution in the DM halo is different 
for different distances from the GC. Today it is commonly be- 
lieved that the profile of a DM halo can be well fitted by the 
universal function 



P;r('") 



(4) 



Here, a is the scale radius of the density profile, which deter- 
mines at what distance from the centre the slope of the pro- 
file changes, e^y and ez are the eccentricities of the DM halo 
within and perpendicular to the Galactic plane and r© is the 
Galactocentric distance of the solar system. In the rest of this 
paper we will assume axisymmetrical haloes, i.e. e^y = 1, but 
is allowed to have values below or above 1, thus allowing oblate 
and prolate haloes. The parameter Pq dm is the DM density of 
the halo at the position of the Sun. The parameters a, /3 and y 
characterize the radial dependence of the density distribution. 
For r a: a the slope of the halo profile is about r for r >> a 
a r"^ and for r << a r Different sets for these parameters 
were suggested by different analyses. Generally the results from 
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Table 1: Parameter settings for the different DM halo pro- 
files considered in this analysis. The parameters of 
the NFW profile are adapted from the publication by 
[Navarro et al. ( 1997) but in this analysis a larger scale 
radius a is used. 



Profile 


a 


/3 


r 


a [kpc] 


Reference 


NFW 


1.0 


3.0 


1.0 


20.0 


Navarro etal. (1997) 


BE 


1.0 


3.0 


0.3 


10.2 


Binnev and Evans (2001) 


Moore 


1.5 


3.0 


1.5 


30.0 


Moore et al. (1999) 


PISO 


2.0 


2.0 


0.0 


5.0 


de Boer et al. (2005) 


240 


2.0 


4.0 


0.0 


4.0 



numerical simulations show that the DM density in the GC is di- 
vergent. Such profiles are called "cuspy" profiles due to the cusp 
in the centre. In this analysis three cuspy p rofiles are considered. 
In the publication bv lNavarro et al.l d 1 9971) it was found that their 
simulation results can be approximated by a profile with a slope 
of y = 1 in the GC, a stronger decreasing slope of /? = 3 for large 
radii and a - 3 for the intermediate distances. This pro file will 
be ref erred to as NFW. From later simulations by Moore et al.l 
d 19991) a profile with a somewhat steeper slope of y = 1 .5 in the 
GC (hereafter Moore) was favoured. A third cuspy halo profile 
was introduced in (Binnev and Evans 2001). They found that the 
microlensing data towards the centre of the MW produced by the 
MACHO collaboration is compatible with a density profile with 
7 > 0.3 (hereafter BE). Also recent N-body simulations showed 
a less cuspy profile for the clumpy component of the DM, which 
is usually called the Einasto profile (Ludlow et al. 2009). The 
amount of matter in clumps is uncertain, but u sually consid- 
ered small compared with the diffuse component (*S pringel et al.l 
l2008 a.b). By considering cusps with slopes between 0.3 and 1.5 
one probably covers the range of possible slopes for cuspy pro- 
files. 

In contrast to cuspy profiles the density distributions pre- 
ferred by observations of rotation curves of low surface bright- 
ness galaxies and dwarf spiral ga laxies have a nearly constant 
DM density in the GC ( y ^ 0) (lOh et al.1 l200l iGentile et al.1 
120071; ISalucci et al.ll2007l) . Such profiles are called "cored" pro- 
files due to the constant density in the kpc scale in the cen- 
tral region. Two different cored halo profiles are considered in 
this analysis. The first profile is called pseudo-isothermal profile 
(hereafter PISO) since it is an isothermal profile (oc r"^) which 
is flattened in the centre. The second cored halo profile (here- 
after 240) is similar to the PISO profile but decreases faster for 
large radii. The parameter settings of the several density pro- 
files are shown in Table [1] The local DM density Pq^dm is a pri- 
ori an unknown parameter and therefore a degree of freedom i n 
the density model . In the publications by lAmsler et al.l (l2008l) : 
iGates et al] (Il995h its value is quoted to be in the range 0.2 - 0.7 
GeV cm"^ (0.005 - 0.018 Mq pc"^). In this analysis Pq^dm is left 
free. In Fig. [T]the different profiles are shown for equal masses 
within the solar orbit. 



3. Dynamical Constraints on the IVIass IVIodel of the 
IVIilky Way 

As mentioned in the introduction the mass model of the Galaxy 
is largely constrained by the total mass of the Galaxy, the rota- 
tion velocity of the Sun at its orbital radius r© and the local matter 
density po.tot - Ptoti^'e), known as the Oort limit. An additional 
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Different halo profiles 




NFW profile (NFW) 




< Binney profile (BE) 




Moore profile (IVIoore) 




PISO profile (PISO) 




240 profile (240) 







Fig. 1: Radial DM density distribution for the parameters given 
in Table [1] The normalization po of the halo profiles is 
fixed for each profile by the requirement of the rotation 
speed of the solar system. 



constraint comes from the local surface density. The experimen- 
tal input for these constraints is discussed first. 

3.1. Total mass 

The total mass of the MW is an important quantity in order to 
constrain the DM density distribution. In general, it is measured 
indirectly either via the kinematics of distant halo tracer stars 
or satellite galaxies or the vertical scale height of the gas dis- 
tribution of the Galactic disc, which can be measured at large 
distances from the GC. 

The definition of the total mass of the Galaxy is difficult 
since a slowly decreasing density has an infinite extension. The 
total mass of a galaxy is conventionally defined as the mass 
within the so-called virial radius. At this radius the total mass 
of the accumulated density of the MW is equal to the mass of a 
homogeneous sphere with the constant density of 200 times the 
critical density of the Universe. 

In the paper by Wil kinson and EvansI d 19991) the mass of the 
MW was estimated from measurements of the radial velocities 
of 27 globular clusters and satellite galaxies for Galactocentric 
distances R > 20 kpc, using a Bayesian likelihood method and a 
spherical halo mass model with a truncated radius. They found a 
mass of the Galaxy within 50 kpc of M(50 kpc) = 5.4+!^ ^ . iqU 

Mq and a total mass of Mtot = 1 -9^; ^ ■ 10^^ Mp. A simi l ar ana - 
lysis was done with more tracer stars bv lSakamoto et alj (l2003h : 
they find Mtot = 2.5;^j'q ■ 10'^ Mq. These measurements used a 
simple parametrization of the potential. Analyses using an NFW 
profile for the DM distribution usually find a lower total mass, 
given it steeper fall-off of the density profile at large distances. 

Using a large sample of 2400 BHB stars from the SDSS in 
the halo (z > 4 kpc, R< 60 kpc) and com paring the results with 
N-body simulations using an NFW profile lXue etaP dlOOSh find 

Mr<6o kpc = 4.0 ± 0.7 ■ 10"Mo, (5) 

which corresponds to Mtot = I-O^qq ' ^O'^ ^p. To get this total 
mass adiabatic contraction (Blum enthal et aljl986h was allowed 
in which case the concentration parameter was taken as a free pa- 
rameter in the fit. This parameter came out to be low (c - 6.6;^['^; 
if adiabatic contraction is ignored, the concentration parameter 
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Comparison of different he 
NFW profile 
BE profile 
Moore profile 

PISO profile 

..... 240 profile 

Luminous matter 

I Critical density 




R [kpc] 



Fig. 2: The mass inside a radius as function of that radius is 
shown for the different halo profiles defined in Table [T] 
The thin solid line represents the visible mass which is 
different for different halo profiles because of the varia- 
tion of the parameters pd and rj. Here the luminous matter 
for the NFW profile is shown. The steep line starting at 40 
kpc represents the mass of a sphere with a density of 200 
times the critical density of the Universe. The crossings 
of the mass distributions with this line represent the virial 
radius and the total Galactic mass of the corresponding 
density distribution. 



becomes around 12 and the mass decreases to Mtot - 0.9 ■ 10'^ 
Mq, which is well within the errors. Similar mass values were 
found from the velocity dispersion by Battaglia et al. (2006), al- 
though with larger errors. These m easurements are co nsistent 
with the favoured ACDM model in i Klvpin et all (|2002|) where 
the virial mass is 1 ■ 10'^ Mq. The value from Eq.|5]will be used 
in the analysis. 

Fig.|2]shows the radial dependence of the total Galactic mass 
for the different spherical halo profiles. At small radii r < 5 kpc 
the density distribution is dominated by the luminous matter, 
shown by the thin solid line, which is independent of the halo 
profile. The mass of a homogeneous sphere with 200 times the 
critical density of the Universe is also shown. The crossing of 
a mass distribution with this line defines the total Galactic mass 
and the virial radius for this density distribution which is 210 
kpc for a Galactic mass of 10'^ Mq. 

The mass distributions of the NFW and the BE profile are quite 
similar since these profiles differ only in the region around the 
GC where the influence of DM is small. The 240 profile yields 
the smallest mass while the mass distribution of the PISO profile 
shows a linear increase with radius. The reason is the quadratic 
decrease (oc r"^) of the PISO profile. Consequently, the integral 
of such a profile leads to a linear increase of the Galactic mass. 
By going from spherical to elliptical profiles the mass can be 
changed significantly, as will be discussed later. 

3.2. Rotation curve 

Each object, which is bound to the MW, is orbiting around the 
GC. Most of the stars and the interstellar medium (gas, dust, 
etc.) are rotating with a velocity distribution v(r). This velocity 



distribution is called the rotation curve (RC) of the MW. 
For a circular rotation of the objects within the Galactic disc the 
rotation velocity is given by the equality of the centripetal and 
the gravitational force 



JO 
dr ' 



(6) 



where v is the rotation velocity at the Galactocentric distance r. 
The gravitational potential O is given by the Poisson equation 



Ad) = -47TGp(r), 



(7) 



where G is the gravitational constant and p(r) is the matter den- 
sity distribution model. For a given DM density profile the gravi- 
tational potential and therefore the rotation curve can be calcu- 
lated and a comparison with experimental data is an important 
check for the DM density profile. 

The rotation curve can be most easily measured from the 
Doppler shifts, like the 21 cm from neutral hydrogen or the rota- 
tional transition lines of carbon monoxide (CO) in the millime- 
ter wave range. A revie w of the various methods was given by 
ISofue and RubinI ( l200lh . A recent summary of all data on the 
rotatio n curve of th e Milk y Way can be found in the publica- 
tion by S ofue et akl (1200 8'). where the numerical values can be 
found on the author's webpage. For radii inside the solar radius 
the distances do not need to be determined, since the maximum 
Doppler shift observed in a given direction is in the tangential 
direction of the circle, so the longitude of that direction deter- 
mines the distance from the center. For radii larger than the solar 
radius the distances to the tracers have to be determined inde- 
pendently, usuall y by die angular thickness of the HI layer, as 
first proposed b v ' Merrifieldl (1 1 9921) . These independent distance 
determinations lead to larger errors in the outer rotation curve. 

The most precise determination and with it the normalization 
of the rotation curve is obtained from the Oort constants, which 
can be determined from the precise distances an d velocities of 
nearby stars, as discussed in most textbooks, e.g. ('B inneviri998l : 
Sparke 2007; Zeilik 1998). These constants are defined as: 



14.4+ 1.2kms"' kpc~' (8) 
-12.0 ±2.8 km s"' kpc"', (9) 



where the experimental values have been taken from 
iKerr and Lynden-Belll (Il986h . One observes that A-B = Vq/iq 
defines the local normalization of the rotation curve, while Ah-B 
defines the slope at the solar position. 

The combination A-B can be more precisely d e termin ed than 
the individual constants. iKerr and Lynden-Belll (Il986l) found 
A-B = Vo/ro = 27+ 2.5 km s"' kpc"'. Using the proper 
mot ion of the black hole in the Galactic centre (Sgr A*) 
Reid and BrunthaleJ (|2004 found 



1 








2 


-ko- 

dr 


re 




1 


dv 






2 


dr 







A-B^ Vo/re = 29.45 ± 0.15 km s"' kpc" 



(10) 



in excellent agreement with recent measurements of paral- 
laxes using the Very Large Baseline Interferometry (VLBI) 
(iReid et all l2009l) . which yield A-B ^We/ie = 30.3 ± 0.9 km 
s~' kpc"'. A further interesting property of the RC is its slope 
at the position of the Sun. This value is strongly connected to 
the values o f A an d B and was determined in the publication by 
iFuchs etalJ ( 120091) to be 



RC Slope 



dln(r) 



A + B 
' A-B 



= -0.006 + 0.016. 



(11) 
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Composition of the rotation curve 




5 10 15 20 25 5 10 15 20 25 

R [kpc] R [kpc] 

(a) Rotation curve ( z = kpc ), non-averaged (b) Rotation curve ( z = kpc ), averaged 



Fig. 3: The rotation curves - calc ulated for differen t halo profiles - in comparison with experimental data, which have been adapted 
from the publication by Sofue et alj ( l2008l) . On the left side all data points are shown, while on the right side a weighted 
average of the experimental data is shown in 17 radial bins. 



From the velocities of stars orbiting Sgr A*, which is con- 
sidered to be the centre of the Galaxy because of its small own 
velocity, the distance between the Sun and the GC has been de- 
termined to (Gillessen et aL,2009,) : 

ro = 8.33 + 0.35 kpc, (12) 

in agreement with previous authors ("Ghez et al.1'2008'). With this 
Galactocentric distance one finds from Eq.[TO]a rotation velocity 
of the Sun 

244+ 10kms"\ (13) 

w hich is consi st ent w ith recent observations of Galactic masers 
in iBovv et aTl (|2009), who used data from the Very Long 
BaseUne Array (VLBA) and the Japanese VLBI Exploration of 
Radio Astronomy (VERA). This speed determines the mass of 
the Galaxy inside the solar radius. 

In this analysis two different rotation velocities are consid- 
ered: the RC within the Galactic disc (Fig. O and the velocity 
distribution for stars outside the disc with z > 4 kpc (Fig. |4|. 
They are discussed separately. 

3.2.1. Rotation curve in tine disc 

For the RC within the disc a combination of different measure- 
ments with different tracers has been summarized by Sofue et al. 
( |2008|) . The experimental data, which can be found on the 
author's web page, were scaled to v© = 244 km s"' at a 



Galactocentric distance of 8.3 kpc. Furthermore, the rotation ve- 
locity was averaged in 17 radial bins from the GC to a radius 
of 22 kpc, as shown in Fig. [3]3 and tabulated in Table |2l The 
shape of the measured velocity distribution shows a strong in- 
crease of the rotation velocity in the inner part of the Galaxy 
which presumably results from the dense core of the Galaxy. For 
the inner Galaxy the rotation curve is dominated by the visible 
matter and the parametrization of section 12.11 yields a reason- 
able description. However, at the outer Galaxy the experimental 
data cannot be explained: all profiles predict a slow decrease of 
the rotation velocity in contrast to the data, which show first a 
decrease between 6 and 10 kpc and then increases again. Such 
a peculiar change of slope cannot be explained by a smoothly 
decreasing DM density profile, but needs substructure, e.g. the 
infall of a dwarf Ga laxy, as mentioned in the introduction and 
Ide Boer et alj (l2005l). Such a ringlike su bstructure is supported 
by the gas flaring dKalberla et al.l 120071) . The thickness of the 
substructure is of the order of 1 kpc, so it should not show up 
for halo stars well above this height; this is indeed the case, as 
shown in Fig.|4] which will be discussed in the next section. 

One may argue about the large uncertainties in the outer rota- 
tion curve, where the distance and the velocity have to be deter- 
mined in contrast to the inner rotation curve, where the tangent 
method yields the distance from the maximum velocity dBinnevI 
Il998h . The rotation curve can be flattened, e.g. by decreasing 
the distance between the Sun and the GC, but then r© ^ 7 kpc 
is needed (Honma and Sofue 1997), which is clearly outside the 
present errors given in Eq.[T2l Also the peculiar change of slope 
near 10 kpc does not disappear. It should be noted that such 
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Circular velocity curve 

NFW profile (10 degrees) 
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(a) Velocity curve (z > 4 kpc), NFW profile 
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(b) Velocity curve (z > 4 kpc) for different DM profiles 
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Fig. 4: The circular velocity curve for halo stars above a height of z > 4 kpc is shown. On the left side the circular velocity curve 
for the NFW profile calculated for different angles with respect to the Galactic disc is shown. On the right side the averaged 
cir cular velocity cur ves for the different halo profiles are shown. The experimental data were obtained from the publication 
bv lXueet all (l2008l) . 



Table 2: Averaged values of the Galactocentric distance and the 
rotation velocity shown in Fig. [3] 



Nr 


fMin 




r [kpc] 


cr, [kpc] 


V [km s ' ] 


(Tv [km s ' ] 


1 


0.0 


1.0 


0.146 


0.002 


281.514 


0.851 


2 


1.0 


1.5 


1.190 


0.011 


261.127 


0.783 


3 


1.5 


2.0 


1.670 


0.016 


247.721 


0.838 


4 


2.0 


2.5 


2.253 


0.016 


236.071 


0.640 


5 


2.5 


3.0 


2.770 


0.020 


233.645 


0.474 


6 


3.0 


4.0 


3.481 


0.017 


242.435 


0.307 


7 


4.0 


5.0 


4.457 


0.021 


257.475 


0.307 


8 


5.0 


6.0 


5.481 


0.026 


257.255 


0.321 


9 


6.0 


7.0 


6.464 


0.027 


261.968 


0.318 


10 


7.0 


8.0 


7.852 


0.001 


251.095 


0.198 


11 


8.0 


8.5 


8.202 


0.009 


243.078 


0.352 


12 


8.5 


9.0 


8.761 


0.046 


233.942 


1.278 


13 


9.0 


10.0 


9.478 


0.052 


231.191 


1.434 


14 


10.0 


11.0 


10.578 


0.199 


238.947 


4.963 


15 


11.0 


12.0 


11.500 


0.227 


253.366 


5.631 


16 


12.0 


16.0 


12.912 


0.115 


253.384 


2.837 


17 


16.0 


22.0 


17.141 


0.796 


274.735 


13.594 



3.2.2. Rotation curve in tine lialo 

The data in Fig.|4]were obtained from a large sample of roughly 
2400 blue horizontal-branch (BHB) tracer stars, as detected in 
Sloan Digital Sky Survey (SDSS), with Galactocentric distances 
up to about 60 kpc and vertical heights of z > 4 kpc. In order 
to connect the observable values - line-of-sight velocity and dis- 
tance - to the circular velocity v(r)= yfFdWfdr the halo star dis- 
tribution function from N-body simulations of the Galaxy with 
an NFW profile was used. 

In order to compare our DM density profiles with these data, 
the velocity curve is calculated for different angles with respect 
to the Galactic disc. Then the results are averaged. In Fig.|4]the 
averaged circular velocity curve and the velocity curves for an 
inclination angle with the normal to the disc of 10°, 45° and 80° 
are shown for the NFW profile. The averaged circular velocity 
curves for the five other spherical halo profiles discussed before 
are shown in Fig.|4j3. The circular velocity distribution is consis- 
tent with the cuspy halo profiles and the PISO profile. The 240 
profile cannot describe the velocity distribution at large radii be- 
cause of the too steep decrease of the density at large radii (oc 



a change of slope happens in other spiral galaxies as well, as 
ca n be seen from the com pilation of rotation curves in Fig. 4 
in lSofue and Rubiiil ( 1200 ll) . For the smooth DM density profiles 
discussed in this paper this feature will be neglected. 



3.3. Surface density and Oort limit 

Jan Oort proposed and performed another interesting measure- 
ment: from the star count as function of their height above 
the disc one obtains the local gravitational potential, which 
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Table 3: Contributions to the local surface density of baryonic 
matter. The total values in the last row include ±lcr 



errors. 



Contribution 


E [Mo pc-^] 


Reference 


Visible stars 


35 ± 5 


Gilmoreet al. (1989) 




27 


Gould etal. (1996) 




30 


Zhena etal. (2001) 


Stellar remnants 


3± 1 


Mera et al. (1998) 


Interstellar gas 


8±5 


Dame (1993) 




13 - 14 


Ollina and Merrifield (2001) 


Total 


35-58 





Vertical potential at the Sun 
NI=W profile 
BE profile 
Moore profile 
PISO profile 
240 profile 




2 2.5 

z [kpc] 



Fig. 5: The vertical gravitational potential at the position of 
the Sun ^(z). The ex perimental data is adapted from 
(iBienavme et al.ll2005h and the errors are assumed to be 
30%. 



is directly proportional to the mass in the plane of the MW. 
Using the precise m e asurem ents from the Hipparcos satellite 
iHobnberg and Flvmil (|2004|) find for the local mass density, 
which includes visible and dark matter. 



P0,tot(z = 0) = 0.102 ± 0.010 Mo pc" 



(14) 



This value was determined by the precise star counts and 
velocity measurements in a volume of 125 pc ar ound the Sun 
by the Hipparcos satellite. iKorchagin et al.l (l2003h analysed the 
vertical potential at slightly larger distances (a vertical cylinder 
of 200 pc radius and an extension of 400 pc out of the Galactic 
plane). For the dynamical estimate of the local volume density 
they obtain the same value with a smaller error: p0,tot(z = 0) = 
0.100 ± 0.005 Mq pc^. To be conservative, we will use the value 
from Eq.[T4] It is called the Oort limit, since it gives the lowest 
value for the local density. 

Integrating the density along the vertical direction within + z 
from the Galactic plane yields the surface density: 



^def (< kl) 



J p(z) dz. 



(15) 



It can be calculated either by integrating the matter density 
distribution directly or by using the gravitational potential of the 



Galaxy. The integration limit is conventionally defined to be 1.1 
kpc. 

According to Eq. ^ this definition is equal to 



^no, (< kl) 



1 

2^ 
1 

2nG 
1 

1 



z 

J A^{r,ip,z') dz' 


dz J \r dr\ dr I dip^ j 



2nG \ dz 



50) 



+ 2(A^ - B^)\z\ 



1 



27rG \ dz 



(16) 



where A and B are the Oort constants discussed before. Since 
they are of the same order of magnitude, the surface density in- 
tegrated to z is proportional to the derivative of the potential at 
height z, as shown by the last approximation in Eq. [16] 

First the surface density of the luminous matter is consid- 
ered. Its experimental value is determined by the summation of 
the different contributions to the luminous matter - the stellar 
population, stellar remnants and the interstell ar gas. A summary 
of the different measurements was given by iNaab and Ostrikerl 
(|2006|) . The surface density of the baryonic matter lies be- 
tween 35 and 58 Mq pc"^ ( Table [3]l, which agrees wit h 48 
+ 9 Mq pc'^ as estimated in iKuiiken and Gilmord (Il991h and 
iHolmberg and Flvnnl (|2004 . 

Unfortunately, our local neighbourhood is not representative 
for the disc, since we live in a local underdensity - the local 
bubble - with an extension of a few hundred pc, which could 
have been caused by a series of rather recent SN explosions 
dMaiz-ApellanidllQOll). Therefore, a fit for the parameters of a 
mass model of the Galaxy might have a somewhat higher surface 
density than the locally observed value. 

The tota l surface density at the posit ion of the Sun was de- 
termined bv lKuiiken and Gilmord (1 19911) to be 



-£{< 1.1 kpc) = 71 ±6Mopc" 



(17) 



from a parametrizati o n of a mass distribution. In the paper by 
Holmber g and Flvnnl (|2004|) the modeling of the vertical grav- 
itational potential resulted in E(< 1.1 kpc) =74 + 6 Mq 
pc~^. The most recent determination of the surface density by 
IBienavme et al.l(l2005h is consistent with both measurements but 
allows somewhat larger errors: 2(< 1.1 kpc) = 68 ± 11 pc"^. 
Comparing this total surface density with the baryonic surface 
density shown in Table [3] it is clear that the local gravitational 
potential is largely determined by the stars, so it will be practi- 
cally independent of the DM halo profile, as shown in Fig. |5]for 
different halo profiles. 



4. Numerical Determination of the IVIass IVIodel of 
the Galaxy 

As mentioned before, the three most important constraints for 
the mass model of the Galaxy are given by the rotation curve 
and the value v©, the total mass Mtot and the local mass density 
Po.tot- This can be easily seen as follows: v| = v^. + v?.|>., where 
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Table 4: Free and fixed parameters for the mass model of the Galaxy and experimental constraints. One observes that there are 5 
free parameters and 8 constraints. Mass densities are in GeV cm"^ or in Mq pc"^, where 1 cm"^ s 37.97 GeV cm"^. 



Free Parameters 



Component 


Parameter 


Symbol 


Value 


Unit 


Halo 


Local DM density 


Pg.dm 




GeV/cm-^ 


Halo 


Scale Parameter 


a 




kpc 


Halo 


eccentricity 






kpc 


Disc 


Density at GC 


Pd 




GeV cm"' 


Disc 


scale length 


rd 




kpc 


Fixed Parameters 


Disc 


scale height 


Zd 


320 


pc 


Bulge 


Density at GC 


Pb 


360 


GeV cm""* 


Bulge 


Eccentricity 


lb 


0.61 




Bulge 


Scale 




0.6 


kpc 


Bulge 


Scale 


ro,b 


0.7 


kpc 


Bulge 


Slope 


Jb 


1.6 




Bulge 


Slope 


Pb 


1.6 




Constraints 


All 


Mass inside 60 kpc 


MR<60kpc 


4.0 ± 0.7 


10" Mo 


Local 


Rotation speed Sun 


Vo 


244± 10 


km s ' 


Local 


Distance Sun-GC 


Tg 


8.33 + 0.35 


kpc 


Local 


Total Surface Density 


S|z!<l.l 


71 ± 6 


Me pc-^ 


Local 


Visible Surface Density 




48 ±9 


Mo pc-^ 


Local 


Mass Density 


Ptot 


0.102 ±0.01 


Mo pc-' 


Local 


Oort Constants 


A-B 


29.45 ± 0.15 


km s"" kpc"' 


Local 


Slope of rotation curve 


d ln(vo)/a ln(r) 


-0.006 ±0.016 





Table 5: Fit results. The units of the different values are given in Table H) The;^^^ contributions are given below the variable value in 
brackets. 







Fitted Parameters 


Derived Quantities 




Nr 


Profile 


PG.DM 


a 


Pi 




CVLI 


Pg.IoI 


Vo 


M60 


M,„, 


A-B 


A(vo) 
dbM 


2l2|<l.l 


Evis 




1 


NFW 


0.32 ± 0.05 


10 


88.3 ± 19.8 


2.5 ± 0.2 
(0.13) 


17.5 


0.094 
(0.58) 


244.5 
(0.0) 


3.9- 10" 
(0.0) 


6.5 • 10" 


29.45 
(0.01) 


-0.002 
(0.06) 


72.2 
(0.04) 


53.8 
(0.41) 


(1.25) 


2 


NFW 


0.27 ± 0.06 


15 


113.1 ± 17.1 


2.4 ±0.1 
(0.01) 


12.0 


0.096 
(0.38) 


244.5 
(0.0) 


4.0- 10" 
(0.0) 


7.4- 10" 


29.45 
(0.0) 


-0.003 
(0.04) 


71.1 
(0.0) 


55.4 
(0.68) 


(1.13) 


3 


NFW 


0.23 ± 0.05 


20 


128.7 ± 69.8 


2.3 ±0.1 
(0.0) 


9.5 


0.097 
(0.27) 


244.4 
(0.0) 


4.1 ■ 10" 
(0.02) 


8.3- 10" 


29.45 
(0.0) 


-0.003 
(0.03) 


70.2 
(0.02) 


56.6 
(0.92) 


(1.27) 


4 


NFW 


0.32 ± 0.04 


10.8 ± 3.4 


91.0 ± 8.0 


2.5 


16.2 


0.095 
(0.50) 


244.4 
(0.0) 


4.0- 10" 
(0.0) 


6.7- 10" 


29.45 
(0.0) 


-0.003 
(0.04) 


72.4 
(0.06) 


54.2 
(0.48) 


(1.07) 


5 


NFW 


0.35 ± 0.06 


14.9 ±4.8 


89.5 ± 8.2 


2.5 


13.4 


0.094 
(0.57) 


244.4 
(0.0) 


4.25 ■ 10" 


9.8- 10" 
(0.0) 


29.45 
(0.0) 


-0.003 
(0.05) 


73.7 
(0.21) 


53.2 

(0.34) 


(1.17) 


6 


NFW 


0.39 ± 0.05 


20.4 ± 4.5 


88.1 ±8.7 


2.5 


11.5 


0.094 
(0.62) 


244.4 
(0.0) 


5.5 ■ 10" 


1.49- 10'- 
(0.0) 


29.44 
(0.08) 


-0.001 
(0.0) 


75.0 
(0.45) 


52.4 
(0.24) 


(1.40) 


7 


NFW 


0.41 ± 0.05 


25.2 ±4.6 


87.3 ± 8.9 


2.5 


10. 1 


0.094 
(0.64) 


244.4 
(0.0) 


6.5 ■ 10" 


1.98- 10'^ 
(0.0) 


29.44 
(0.0) 


-0.001 
(0.09) 


75.9 
(0.68) 


51.9 
(0.19) 


(1.61) 


8 


BE 


0.25 ± 0.05 


10.2 


133.6 ± 15.2 


2.29 ± 0.09 
(0.0) 


17.6 


0.096 
(0.31) 


244.4 
(0.0) 


4.1 ■ 10" 

(0.01) 


7.5 • 10" 


29.45 
(0.0) 


-0.003 
(0.03) 


70.6 
(0.01) 


56.2 
(0.84) 


(1.20) 


9 


Moore 


0.25 ± 0.05 


30.0 


114.7 ± 17.3 


2.36 ±0.11 
(0.01) 


6.2 


0.095 
(0.32) 


244.4 
(0.0) 


4.1 ■ 10" 
(0.01) 


7.6- 10" 


29.45 
(0.0) 


-0.003 
(0.04) 


70.4 
(0.01) 


56.2 
(0.82) 


(1.21) 


10 


PISO 


0.20 ± 0.04 


5.0 


150.4 ± 12.9 


2.21 ±0.08 
(0.02) 


46 


0.098 
(0.19) 


244.4 
(0.0) 


4.1 ■ 10" 

(0.03) 


1.45 • 10'^ 


29.45 
(0.0) 


-0.004 
(0.02) 


69.3 
(0.08) 


57.8 
(1.18) 


(1.54) 


11 


240 


0.26 ± 0.03 


4.0 


53.1 ±9.6 


3.0 ±0.2 
(1.36) 


26.3 


0.095 
(0.51) 


244.6 
(0.0) 


1.7 ■ 10" 

(10.45) 


1.8- 10" 


29.47 
(0.02) 


-0.002 
(0.06) 


70.0 
(0.03) 


55.0 
(0.61) 


(13.04) 


12 


NFW 


0.52 ± 0.07 


26.6 ± 4.9 


84.2 ± 9.5 


2.5 


9.6 


0.094 
(0.66) 


244.3 
(0.01) 


6.6- 10" 


1.95 • 10'' 
(0.03) 


29.43 
(0.02) 


-0.001 
(0.10) 


80.3 
(2.38) 


50.1 
(0.66) 


(3.24) 



v^j^ and Vpj^ are proportional to pvis and poM, respectively; for 
a given halo profile Mtot is determined by pdm, while the Oort 
limit po,tot determines po.tot = Pvis + Pdm- So in principle one has 
3 constraints with only 2 variables pvis and pdm, //the shapes of 
the DM halo and the visible matter would be known. 

Unfortunately, additional important parameters are i) the ec- 
centricity of the DM halo ii) the concentration of the DM halo 
iii) the scale length of the disc iv) the mass in the bar/bulge. 



In addition the mass model is sensitive to the geometry, i.e. the 
Galactocentric distance from the Sun d© and the halo profile. 
Additional constraints come from the surface density, but here 
the visible surface density has a large uncertainty as discussed 
before. The parameters and constraints have been summarized 
in Table |4] The parametrization of the mass of the bulge was 
chosen to describe the rotation curve at small radii, which works 
reasonably well, as can be seen from Fig. [3] Given that the mass 
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Fig. 6: The local DM densities Pq.dm are shown for different fits with different parameters. The numbers correspond to the numbers 
of the fit results in Table |5] 



model is not very sensitive to this inner region, the parameters 
of the bulge will not be varied anymore. 

To optimize the remaining parameters in order to best de- 
scribe the data, the following function w as minimized using 
the Minuit package dJames and Roos|[l975l) 

2 ' 2 ' 2 

rr (T cr 

(^Qt.-^)' ((A-Br'^-D)2 

2 ' 2 ^ 

The index means the observables were calculated from the 
fitted parameters, while D denotes the experimental data for the 
observable and cr its error. The constraints have been summa- 
rized in Table |4] 

The fit shows a more than 95% positive correlation between 
the local dark matter density and the scale length of DM halo a 
and an equally large negative correlation with the scale length rj 
of the baryonic disc. Consequently, it is difficult to leave param- 
eters free in the fit. Therefore the fit was first performed for fixed 
values of a (rows 1-3 of Table |5]l and then rd was fixed (rows 4- 
7). With the other free parameters all experimental constraints 
could be met, as indicated by the;^^ values in brackets below the 
fitted values in Table|5] Of course, the total mass changed for the 
different fits. Fig.|6]shows the resulting local DM density versus 
the total mass, as calculated from the fitted parameters. It shows 
that in spite of the small errors for the local density in individual 
fits the spread in density is still quite large. 

The fit was repeated for other halo profiles, which gave simi- 
larly good;^'^ values, as shown by rows 9-11 in Table |5] So with 
the present data one cannot distinguish the different halo profiles. 

Sofar only spherical halos have been discussed. Allowing 
oblate halos with a ratio of short-to-long axis of 0.7 the local 



DM density increases by about 20%, as shown by the last row 
of Table |5] As mentioned before, dark discs can enhance this 
value considerably more, so the uncertainty usually quoted for 
the local dark matter densit y in the range of . 2 to 0.7 GeV cm"^ 
(0.005 - 0.018 Ma pc-3) (lAmsler et alJl2008t iGates etanfT99l 
is still valid in spite of the considerably improved data. 

5. Conclusion 

In this analysis five different halo profiles are compared with 
recent dynamical constraints as summarized in Table |4] The 
change of slope in the RC around 10 kpc (Fig.O was ignored, so 
the monotonical decreasing RC for the smooth halo profiles do 
not describe the data well. The change of slope may be related 
to a ringlike D M substructure , as indicated by the structure in 
the gas flaring (Kalberla et al."2007) and by the structure in the 
diffuse gamma radiation (de Boeret al. 2005). Such a ringlike 
structure of DM gives a perfect description of the rotation curve, 
especially the fast decrease between 6 and 10 kpc. If the DM sub- 
structure is included, the local DM density increases above the 
values found in this analysis, so the values quoted here should 
be considered lower limits. 

The astronomical constraints are consistent with a density 
model of the Galaxy consisting of a central bulge, a disc and an 
extended DM halo with a cuspy density profile and a local DM 
density between 0.2 GeV cm^^^ (0.005 Mq pc^^^) and 0.4 GeV 
cm"^ (0.01 Mq pc^^*), as shown in Fig. |6] Strong positive and 
negative correlations between the parameters were found in the 
fit and they are causing the obvious correlations between Pq.dm 
and Mtot in Fig.|6] For non-spherical haloes these values can be 
enhanced by 20%. If dark discs are considered, densities up to 
0.7 GeV cm"-' (0.018 Mq pc^^) can be easily imagined, so the 
previous quoted range of 0.2 - 0.7 GeV cm^^ (0.005 - 0.018 Mo 
pc"^) seems still valid. This range is considerably larger than the 
values quoted by analyses which used a Markov Chain method 
to min imize the likelihood; th ey find Pq.dm = 0.39 ± 0.03 GeV 
cm-3 (ICatenaand UUiol l2009l) and po.dm = 0.32 + 0.07 GeV 
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cm~^ (IStrigari and Trottall2009h respectively. But given the good 
values for our fits obtained for a large range of DM densities 
we see no way that the errors can be as small as quoted by these 
authors. 
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